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Using a chirally invariant effective Lagrangian, we calculate the density and isospin dependences of 
the in- medium axial coupling, g A , in spatially uniform matter present in core collapse supernovae 
and neutron stars. The quenching of g\ with density in matter with different proton fractions is 
found to be similar. However, our results suggest that the quenching of the nucleon's g* A in matter 
with hyperons is likely to be significantly greater than in matter with nucleons only. 
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The accurately measured beta decay lifetime of the 
neutron in vacuum, n — ► p + eT + D e , fixes the ratio 
of the axial and vector couplings of the neutron to be 
\g A /g v \ = 1.2601 ± 0.0025 §. Studies of beta decays 
in nuclei, however, have long suggested that a value of 
\gA/gv\ — 1 better fits the observed systematics [Q. 
Such a lower value also appears to be consistent with 
pion-nucleus optical potentials and the systematics 
of Gamow- Teller resonances in nuclei (i). Data from 
muon capture on nuclei, in which the relevant momen- 
tum transfer q 2 ~ — 0.9m^, have been recently analyzed 
including detailed nuclear structure effects || with the 
conclusion that a quenched gA (assuming gy = 1) is not 
necessary inasmuch as the vacuum value of gA adequately 
accounts for the data. 

The above experiments measure space- like axial transi- 
tions in nuclei. At finite density, however, space-like and 
time-like axial matrix elements are not necessarily equal, 
since Lorentz invariance is broken. In fact, there are in- 
dications from experiments with first-forbidden (3 decays 
of light nuclei that the time-like axial charge increases 
by about 25% in medium g. A theoretical expectation 
of this enhancement in terms of soft-pion exchanges has 
been offered in Ref. J7J. 

The space-like quenching of gA in nuclei at low mo- 
mentum transfers has been attributed to a combination 
of effects, including the partial restoration of chiral sym- 
metry in a nuclear medium, the direct participation of 
the A(1232) in renormalizing the in-medium axial- vector 
current, and tensor correlations in nuclei in which shell 
structure effects are important An illuminating dis- 
cussion of the extent to which the quenching phenomenon 
is intrinsic to the basic property of the "vacuum" de- 
fined by a baryon-rich medium has been given by Rho 
(cf. Ref. H and references therein). Later discussions of 
the quenching phenomenon in chiral approaches to spa- 
tially uniform matter can be found in Ref. |5). The is- 
sue of breaking and restoring fundamental symmetries 
at large baryon density is presently intractable in lat- 
tice gauge simulations. We therefore employ an effective 
field-theoretical approach, based upon chiral symmetry, 
to consider medium modifications of the nucleon's axial 
coupling. 



The precise value of the in-meidum axial coupling, de- 
noted by g* A hereafter, in the dense matter encountered 
in astrophysical phenomena such as core collapse super- 
novae and neutron stars is crucially important. In these 
cases, weak interaction rates (that are oc g* A 2 ) drive the 
dynamics from beginning to end. In contrast to labora- 
tory nuclei, in which the weak processes occur at nuclear 
to subnuclear densities, astrophysical settings feature su- 
pernuclear densities. This highlights the need for knowl- 
edge of g* A well beyond nuclear densities and in uniform 
matter with varying isospin content. In the supernova 
environment, e _ -capture reactions on neutrons and nu- 
clei begin the process of neutronization and decrease of 
the total lepton fraction, Yl — n^/riB, whose value af- 
ter ^-trapping (~ 0.38 — 0.4) determines the masses of 
the homologous core and initial proto-neutron star and 
thus the available energy for the shock and subsequent v- 
emission [jl~o|| . The bremsstrahlung (n+n — > n+n + v+v) 
and modified Urea (n + p — > n + n + e + + v + v) processes 
dominate in the production and thermalization of /z and 
r neutrinos |TjJ . The (/-luminosity and the time scale 
over which vs remain observable from a proto-neutron 
star are also governed by charged and neutral current 
interactions involving baryons at high density fl^| . The 
long-term cooling of a neutron star, up to a million years 
of age, is controlled by i/-emissivitics from the densest 
parts of the star; thereafter the star is observable through 
photon emissions, which may allow us to determine the 
star's mass, radius, and internal constitution ]l3] ]. 

In this Letter, we focus on the behavior of the in- 
medium space-like g* A , both with increasing density and 
with varying isospin content, in spatially uniform mat- 
ter in which shell effects characteristic of finite nuclei 
are absent. This enables us to expose cleanly the role 
played by chiral symmetry alone. Our results are di- 
rectly relevant to astrophysical situations in which spa- 
tially homogeneous matter are encountered. We consider 
the cases of isospin symmetric matter (proton fraction 
x = n.p/riB = 1/2), pure neutron matter (x = 0), and 
neutron-star matter in which the equilibrium value of 
the proton fraction, i, is determined by the conditions 
of beta stability and charge neutrality. The latter two 
cases have been largely ignored in the literature, but are 
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essential for astrophysical modeling in which neutrinos 
play a dominant role. 

The nucleon coupling to the axial current is directly 
determined through the matrix element 

(N(p 2 )\A p \N( Pl )) =u(p2)^T[ ltll5 g A 

+q^i5hA(q 2 )] u{jpt) , (l) 

where the u(j>2) and u(pi) are the spinor solutions of 
the Dirac equation for nucleons and q = pi — Pi ■ In this 
form, gA is the axial coupling and hA{q 2 ) the form factor. 
Taking medium effects into account does not change this 
definition when we consider a medium- modified g* A . In 
the applications of interest here, momentum and energy 
transfers are moderately low. We therefore concentrate 
on medium modifications of gA [@ • The axial current 
can be computed from first principles. Explicitly, 
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where the sum runs over all fields in the Lagrangian and 
<5g<I>i is the change in each field under an axial trans- 
formation. The fact that the axial coupling is related 
to pion dynamics requires that our approach should in- 
volve the spontaneous breaking of chiral symmetry. We 
therefore use a chiral Lagrangian to describe the inter- 
actions between baryons, keeping in mind that it should 
reproduce the vacuum value of gA, equilibirum empirical 
properties of isospin symmetric and asymmetric matter, 
reasonably describe closed-shell nuclei, and account for 
observables in pion-nucleon interactions. Fortunately, it 
is possible to achieve these objectives with mean fields 
computed at the tree level. Originally motivated by the 
gluonic trace anomaly in QCD, a model incorporating a 
heavy glueball field was devised which replaced the chiral 
"sombrero" potential with a logarithmic one Jl5|,[l6| . For 
our purposes, it will be sufficient to "freeze" the glue- 
ball field at its vacuum value, which shall have negligible 
consequences since the gluon condensate changes only 
slightly at finite density |L7||. Specifically, we have 
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The field strength tensors are 



F, 



d p p v - d v p p + g p p p x p„ + g p a p x a v , 
G pv = d p a v - d u a p + g p p p Xa„- g p p u x a p , 



with the covariant derivatives 

A M cr = d p u + g p a p -ir , 

A p ir = d p 7? + g p p p x if - g p aa p . 

The N are isospinor nucleons, a and 7r are scalar and 
pseudoscalar isovector chiral mesons, and p and a\ are 
isovector mesons which are the vector and psuedovector 
representations of a gauged chiral symmetry group and 
make vector meson dominance inherent in the model p8[ . 
The a; is a chiral singlet vector field which supplies the 
short-ranged repulsion necessary for saturation, with a 
self-interaction term included to soften nuclear matter. 
The non-standard coupling between the nucleon and me- 
son fields, with the coefficient D, serves to generate the 
physical value of gA at the mean-field level while respect- 
ing chiral symmetry ||l9|| . Note that the pion-nucleon in- 
teractions in Eq. (m are chiral invariant. 

The effective chiral Lagrangian of Eq. (||), like that 
of the Walecka model and its variant approaches to nu- 
clear matter ^0|, is to be used at the mean-field level. 
This model is able to describe closed-shell nuclei similar 
in accuracy to that obtained in Quantum Hadrodynami- 
cal models [ 17j . Calculations of low-energy pion-nucleon 
scattering also compare well with data . The scaling of 
the physical constants, effective masses, and the quench- 
ing ofg^ in isospin symmetric matter are also detailed in 
Ref. ||. Note also that the Goldberger-Treiman relation, 
9-nNNf-n = 9aM, remains valid at finite baryon density, 
since the model encodes chiral symmetry. 

We turn now to spatially uniform matter with an 
unequal number of neutrons (n) and protons (p). In 
medium, the a, ujq, and pfj fields develop non-zero ex- 
pectation values, which we denote by <7, u>, and p. These 
are obtained by solving the equations of motion: 
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= g(NN) =g{n; + <) , 



mlto + 4G 4 w 3 = g u (rtN) = g u (n p + n„) , 
TO pP = 9 p (N^t 3 N) = g p (n p - n n ) 



(4) 



Above, n\ and m with i = n,p denote the scalar and 
baryon number densities, and r 3 is the third component 
of the isospin operator. 
From the Dirac equation, 



ga) N = 0, 



(5) 



we find the effective mass, M* = ga, and the effective 
chemical potentials, p* = pi — g^oJo — ^9 p t 3 Pq- The 
scalar and baryon densities depend self-consistently on 
these effective parameters: 
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where E* f . = ^k 2 h + M* 2 and k fi = \/ li* 2 - M* 2 . In 

pure neutron matter, n p = and = 0. In beta-stable 
and charge neutral neutron-star matter, the constraints 

n p - n e - n M = , (7) 

where n e and n M are the number densities of electrons and 
muons, determine the equilibrium proton fraction. With 
the solutions to Eqs. (|]) the energy density of matter, 

£ = \mlu? + G 4 tl> 4 + \m 2 p p 2 + \B - In ^ - l) 

+ J2^J h dkk 2 y/k 2 +M* 2 , (8) 

and the energy per baryon, E/A — £ /ns — M, are easily 
computed. 

In order to compute the axial current, we will need 
only the terms in Eq. ^ which are finite for a nonzero 
sigma mean field and contribute to the nucleon matrix 
element Eq. ([[]). The relevant terms in the Lagrangian, 
Eq. (^|), will be the nucleon kinetic energy and deriva- 
tive coupling to the pion, the latter since S^TT b = S ab cr. 
Renormalization of n — a\ mixing introduces the physical 
fields 

*-{{-%)*> 1- w 

where the effective mass of the a\ meson depends on the 
sigma mean field as m* 2 = m 2 + (m 2 — m 2 )<7 2 /a 2 . This 
replacement leads to additional derivative terms from the 
ai-nucleon coupling and, in terms of the redefined pion, 
contributions to the matrix element arise from 



with its mean-field dependence made explicit. Density 
dependence is thus generated by the scalar mean field, 
a, alone. We note that the functional form of Eq. (|T^), 
orginally derived in Ref. j^] for isospin symmetric mat- 
ter, remains intact for isospin asymmetric matter. The 
numerical value of g A , however, is controlled by the mag- 
nitude of a which depends on the proton fraction of mat- 
ter. This relation and its implications to be discussed 
below are among the prinicpal results of this work. 

We now specify the physical parameters in Eq. (|l|) 
and discuss quantitative results in matter. Fixing the 
vacuum value of the nucleon mass determines the chiral 
coupling g = 9.2. We take ct = U/VZ^ = 102 MeV, 
where is the pion renormalization constant necessi- 
tated by 7r — a\ mixing. The value D — 1.17 reproduces 
gA = 1.26 in vacuum. The saturation of nuclear mat- 
ter at its empirical density jiq = 0.15 fm~ 3 and energy 
per baryon E/A = -16 MeV require B = (323 MeV) 4 
and g u = 11.4. The omega self-interaction strength is 
G = 0.19g w , fixed to produce maximal softening of the 
EOS while not generating spurious field solutions. At 
saturation, the resulting effective mass is M* = 0.68M 
and the compression modulus K = 320 MeV. The rho 
coupling g p = 8.0 yields the empirical symmetry energy 
of =30 MeV. 

The procedure adopted above to fix the various cou- 
plings in Eq. ([j]) implies that the behavior of g* A with 
increasing density, including its value at the equilibrium 
density no of symmetric nuclear matter, is to be regarded 
as a prediction, albeit within the confines of the model 
adopted for matter. The EOS of neutron-star matter at 
supra-nuclear densities is subject to the constraint that 
it must support at least 1.44M , which is the most accu- 
rately measured mass of the neutron star in the binary 
pulsar PSR 1913+16 pl| . Presently, more severe con- 
straints at high density are not available. 
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where the medium-dependent renormalization constant 
Z* = l- g 2 a 2 /m* a 2 . Using Eq. (|), we find 
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Inserting this into Eq. ([!]) we have, in terms of the sigma 
mean field and vacuum constants, 



g A =\l + D- 



9> 2 



m 2 + (m 2 — m 2 /) a 2 /a 2 



(12) 



Above, the first factor contains the "bare" axial coupling, 
including the standard axial interaction and the mean- 
field dependent modification from the D term. The sec- 
ond factor is the pion renormalization constant of Eq. (j^) 
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FIG. 1. Left panel: The energy per baryon in nuclear 
(dotted lines), pure neutron (dashed) and beta-stable neu- 
tron-star matter (solid) as functions of baryon density in units 
of no = 0.15 fin -3 . Right panel: Particle concentrations in 
beta-stable neutron star matter. 
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The left panel of Fig. [l] shows how E/A varies with 
ns/no as the proton fraction x is varied from nuclear to 
beta-stable neutron-star to the ideal case of pure neutron- 
star matter. The right panel shows the concentrations 
Yi = jii/riB of n, p, e~, and fjT in neutron-star mat- 
ter. The corresponding EOS yields a maximum mass of 
M max /M.Q = 1.6; for this configuration the central den- 
sity n c /no = 7.4 and radius R — 10.24 km. These values 
are to be compared with M mal /M = 2.2, n c /no = 5.3, 
and R = 12.74 km for the case of pure neutron matter, 
which highlights the role of the isospin content in matter 
for this EOS. 

Fig. H shows the density and proton fraction depen- 
dences of the scalar mean field a/ao (left panel) and the 
axial- vector coupling g* A (right panel) . The quenching of 
both of these quantities with increasing ub is clearly ev- 
ident, with a depreciation of 12% at no and a 19% drop 
at 4no- The relatively mild variation with x is important 
insofar as guidance for the quenching of g A from labora- 
tory studies ( which sample a narrow range in x) of nu- 
clei have the potential of being directly and immediately 
useful in astrophysical applications (in which a broader 
range of x is sampled). For modeling purposes, a simple 
parametrization of g* A in terms of baryon density is 

g* A ~g A (l-—£*—J). (13) 
V 4.15 (n + n B )/ 

This expression matches the results of Eq. ( p"2"| ) to within 
1% accuracy for all Ub < 4.5no- 
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FIG. 2. The sigma mean fields (left panel) and the axial 
coupling constants (right panel) in nuclear (dotted lines), pure 
neutron (dashed), and beta-stable neutron-star matter (solid) 
as functions of baryon density in units of no = 0.15 mi . 

The quenching of gA considered in this work stems 
chiefly from the medium-dependent scalar field and ir—ai 
mixing. The combination of these two effects has been 
shown to be virtually independent of the isospin con- 
tent, suggesting that the in-medium behavior observed 
in nearly iso-symmetric matter will be present in neu- 
tron and stellar matter as well. The excitation of the A, 
which is crucial in non-relativistic descriptions of quench- 



ing, has yet to be satisfactorily implemented in a rel- 
ativistic field theory. Its addition would likely lead to 
further reduction. 

It is also worthwhile to point out here that a falls more 
rapidly with density in matter with hyperons than with- 
out hyperons (see Fig. 6 of Ref . j2^] ) . This is mainly due 
to the presence of additional baryonic components with 
dissimilar masses. Consequently, the nucleon's g* A would 
be quenched to a greater extent in the presence of hyper- 
ons. A verification of this expectation would require an 
extension of flavor symmetry, which has been attempted 
recently with only limited success p3[ . Our results in 
this work suggest that the extension of the chiral model 
to incorporate strange and Delta resonances with the full 
effects of relativity would be worthwhile. 
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